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DE-EMBEDDING  MILLIMETER-WAVE  INTEGRATED  CIRCUITS  WITH  TRL 


I.  INTRODUCTION 

Now  that  you  have  the  use  of  a  network  analyzer  that  goes  to  40 
GHz  or  above,  just  what  do  you  do  with  it?  Using  the  calibration 
procedures  included  with  the  analyzer  will  let  you  look  at  pack¬ 
aged  millimeter-wave  circuits.  The  two-port  parameters  will 
include  the  effects  of  line  and  transition  reflections  and  loss¬ 
es.  Using  the  line  extender  artifice  is  useful  only  if  the 
transitions  are  near  ideal  and  uhe  line  parameters  well  known. 

De-embedding  to  the  circuit  or  device  contained  within  the  line 
(microstrip,  coplanar,  suspended  substrate,  etc.)  usually  re¬ 
quires  ideal  or  measured  shorts,  opens,  and  resistive  loads.  At 
normal  microwave  frequencies  a  short  may  appear  as  such,  but  at 
millimeter  wavelengths  a  short  will  surely  be  inductive,  an  open 
radiative,  and  a  resistor  inductive.  At  millimeter  wavelengths 
one  very  interesting  de-embedding  procedure  called  TRL1 
(Through-Reflect-Length)  seems  almost  magic-like  in  that  very 
non-ideal  elements  can  be  used  for  the  procedure.  These  elements 
are: 

T:  Through  line  which  connects  the  two  transitions. 

R:  Reflect  elements  which  can  be  just  about  anything  as  long  as 

identical  ones  are  used  to  terminate  the  input  and  output 
transitions . 

L:  Length  of  transmission  line  connecting  the  two  transitions. 


TRL  de-embedding  is  easy  to  implement  on  wafers  which  can  use 
coplanar  probes2  or  microstrip  launchers3'4.  The  input  and 
output  probes  can  be  moved  apart  thereby  facilitating  the  "L" 
part  of  TRL. 

The  method  to  be  described  will  allow  TRL  de-embedding  using  any 
type  of  transmission  line  without  the  use  of  external  probes.  The 
abilty  to  use  waveguide  input  and  output  transitions  should  be  of 
particular  interest  to  designers  of  millimeter-wave  hybrid  inte¬ 
grated  circuits. 


II.  METHOD 

Figure  1  indicates  the  four  measurements  necessary  for  de-embed¬ 
ding.  Note  that  TRL  transition  A  and  TRL  transition  B  both 
include  line  lengths  which  meet  somewhere  in-between  (reference 
position) .  The  actual  transitions  which  connect  to  the  analyzer 
ports  can  be  quite  different,  and  may  even  include  difierent  bias 
lines . 
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Figure  la.  shows  the  basic  through  structure.  S  parameters  are 
taken  over  the  desired  frequency  range.  The  transitions  can 
actually  have  poor  VSWR  without  upsetting  the  de-embedding  proce¬ 
dure. 

Figure  lb.  shows  the  reflect  measurement  used.  Since  any  two 
identical  reflections  may  be  used,  it  is  particularly  simple  to 
remove  equal  large  sections  from  both  sides  of  the  reference 
position.  The  reflection  elements  are  therefore  reduced  sections 
of  line  terminated  with  a  radiative  open  circuit. 

Figure  lc.  shows  the  length  measurement  configuration.  Note  that 
TRL  transition  A  and  TRL  transition  B  are  the  same  as  in  the 
other  measurements  except  for  an  added  length  of  transmission 
line.  This  requires  another  fixture  exactly  like  that  of  Fig. 
la.  except  for  an  extended  length  in  the  region  surrounding  the 
reference  position.  To  insure  "exactness”  the  fixtures  are  best 
made  on  a  computer-controlled  mill. 

Figure  Id.  is  the  same  as  figure  la  except  that  the  device  to  be 
measured  is  shown  centered  at  the  reference  position. 

The  four  sets  of  measurement  data  should  be  stored  in  data  files 
for  manipulation  by  the  de-embedding  program.  The  mathematics 
behind  the  procedure,  based  (with  exceptions)  on  Reference  1,  is 
given  in  Appendix  A.  A  de-erabedding  program  which  directly  fol¬ 
lows  this  is  given  in  Appendix  B. 
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Fig.  la  I - 1 
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Fig.  lb 
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Fig.  Id 


D  U  T 


DEVICE  MEASUREMENT 


Figure  1.  The  four  measurements  required  for  TRL  de-embedding 
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III.  An  Example  of  De-Embedding 

To  illustrate  the  effectiveness  of  the  above  procedure,  two  Ka 
band  suspended  substrate  fixtures  with  waveguide  transitions  were 
fabricated  using  a  computer-controlled  mill5.  The  fixtures  were 
fabricated  for  line  lengths  (excluding  probes)  of  1.4  and  2.0 
inches.  The  length  difference  is  about  twice  the  suspended 
substrate  wavelength  at  33  GHz.  Instead  of  using  .080"  waveguide 
probe  lengths  which  would  have  given  excellent  VSWR  for  the  92 
ohm  transmission  lines6,  .050"  probes  were  used.  A  typical 
suspended  substrate  bandpass  filter  was  inserted  in  the  line  as 
the  DUT.  Figures  2  and  3  show  the  detrimental  effects  of  using 
the  intentionally  bad  probes. 

The  two-port  S  parameters  for  the  above  two  cases  were  measured 
in  the  smaller  fixture  as  well  as  the  S  parameters  of  a  non-ideal 
open  circuit  (about  l  inch  of  line  removed  in  the  center) .  Next 
the  S  parameters  of  the  2 -inch  line  were  measured  in  the  larger 
fixture.  All  of  the  measured  parameters  were  put  into  computer 
files  which  were  then  called  up  by  the  de-embedding  program.  The 
results  of  de-embedding  the  filter  are  shown  in  Figure  4.  The 
usual  flat  top  and  reasonable  return  loss  are  similar  to  measure¬ 
ments  made  with  normal  .080"  probes. 


IV.  Conclusions 

The  use  of  two  circuit  fixtures,  one  larger  than  the  other  by  a 
nominal  length,  leads  to  a  de-embedding  procedure  which  should  be 
very  useful  at  millimeter  wavelengths.  Particularly  important, 
any  type  of  transmission  line  may  be  used.  Transitions  and  bias 
lines  may  be  different  for  the  input  and  output.  A  transistor 
chip  with  bond  wires  may  be  measured  in  the  same  media  as  its 
intended  use.  Discontinuities  are  therefore  measured  along  with 
the  chip  and  can  be  used  directly  in  subsequent  circuit  design. 
The  reproducibility  of  the  fixtures  and  transmission  lines  will 
determine  the  upper  frequency  limit.  Using  .0005-in.  tolerances, 
useful  parameters  should  be  obtainable  to  at  least  100  GHz. 


4 


GHz  -*•  GHz 


j  (jure  3.  Measured  parameters  of  bandpass  filter  figure  4.  be-embedded  filler 

with  poor  transitions 


Appendix  I:  TRL  De-Embedding 
A.  FINDING  THE  LINE  PARAMETERS 

Both  the  attenuation  and  effective  dielectric  properties  of  the 
transmission  line  may  be  found  from  the  measurements  indicated  in 
figures  la  and  lc.  If  there  were  no  transition  reflections,  then 
it  would  be  simple  to  get  these  values  from  the  two  S21  measure¬ 
ments.  However,  reflections  from  the  transitions  introduce 
errors  which  can  be  eliminated  by  the  first  part  of  a  TRL  pro¬ 
gram. 

The  scattering  parameters  derived  from  each  of  the  four  steps  can 
easily  be  converted  to  transmission  parameters: 


Scattering  Parameters 

bi=  s11a1  +  s12a2 
b2  =  s21al  +  s22a2 

The  relationship  between  the 
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where  ds=Slls22  -  s12 


Transmission  Parameters 

bi:rlla2  +  r12b2 
al-r21a2  +  r22b2 


two  parameters  is 
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and  the  "R"  matrix  used  here  is  sometimes  called  the  "transmis¬ 
sion"  matrix. 


The  R  matrix  has  a  property  similar  to  ABCD  matrices;  they  may  be 
multiplied  in  sequence  (cascaded) . 
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Measure 


1 . 


SAB 


the  S  parameters  of  the  through 
connection. 
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Convert  to 


AB 


noting  that 


RAB 
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2 .  Measure 


SALB 


the  S  parameters  of  the  through 
connection  with  the  added  length 
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ralb 
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= 

ra 
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There  is  enough  information  from  these  two  measurements  to  obtain 
the  parameters  of  the  transmission  line  used,  i.e.,  the  attenua¬ 
tion  and  phase  velocity  at  each  frequency.  The  added  length  must 
have  the  same  characteristics  as  the  line  in  which  the  DUT  is  to 
be  embedded. 


3.  Find  the  matrix 


[ralb 


(1) 


Note  that  all  the  elements  of  the  T  matrix  are  known  since  the 
elements  of  RAtB  and  R^yg  were  directly  found  from  the  SALB  and 
SAB  measurements.  *  * 

Using  the  cascade  properties  of  the  R  matrix  we  have 
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so  that 


T 


-1 


(2) 


The  added  line  of  length  L  is  non-reflecting  (same  impedance  as 
end  of  TRL  transitions) ,  so 


where  r=  a  +  j/3 

a=  attenuation  const 
)S=  propagation  const 


Tneref ore 


Using 


we  have 


T1  T2 
t3  t4 


RA1  RA2 
RA3  RA4 


RA1  RA2 
RA3  RA4 


-rL 


0 

,rL 


Multiplying  matrices  we  find: 


(3)  T1.RA1  +  T2.RA3 

-  RAle'TL 

(5) 

TI ‘ RA2 

+  T2 ' RA4  ~  RA2eTL 

(4)  T3.RA1  +  T4.RA3 

=  KA3e  TL 

(6) 

T3  *  RA2 

+  T4.RA4  =  RA4eTL 

From  3,4  we  have 

rai 

T2 

e'TL  -  T4 

RA3 

e_ 

T  L  m 

T3 
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so  that 


0 


(8) 


The  solutions  for  e  TL  and  erL  are  the  two  solutions  to  the 
complex  equation 

G2  -  G(T1+T4)  +  (TxT4  -  T2T3)  =  0 


we  find 

G  =  B*  j^l  ±  D1/2  j  with  solutions  |G1|*e3tgl,  | G2 1  * 

Tx  +  T4  C 

where  B  =  -  C  =  T1T4~T2T3  D  =  1  -  - 

2  B2 

The  solution  for  e-TL  =  e_aL* e”j^L  will  be 

Ig-jJ  if  tgl  is  negative 

|G2|.e^g2  if  tg2  is  negative 

The  attenuation  (dB)  in  the  distance  L  is  20*log10(G1  2)  so 

(12) 

where  G3  ,  2  is  the  correct  solution  as  found  above 
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27r*L*f  *yeeff 

Also,  since  /3L  =  tg  =2ir  •  L/wavelength  =  - 


eeff  =  {  30  •  tg/ [2 . 54  •  2jt  •  L(in)  -f  (gHz]  }2 


(13) 


We  have  found  the  transmission  line  parameters  dB/inch  and  eeff 
from  measuring  the  S  parameters  of  the  two  fixtures  with  differ¬ 
ent  lengths. 

For  low  loss  media  such  as  suspended  substrate,  the  added  length 
L  should  be  large  enough  so  that  the  transmission  loss  of  the 
larger  line  will  be  evident.  The  added  line  may  need  to  be  a 
centimeter  or  longer.  Since  measurements  of  the  phase  of  sll  and 
s22  are  ambiguous  by  multiples  of  2 v  radians  if  L  is  larger  than 
1/2  wavelength,  one  can  resolve  the  ambiguity  in  tg  as  follows: 

a.  Estimate  eef£  and  find  an  approximate  tg. 

b.  Use  this  value  to  calculate  eef^°,  eefj+,  and  eeff“  for  the 
approximated  values  of  tg,  tg+27r,  ana  tg-27r  respectively. 

c.  Choose  which  of  the  three  values  comes  closest  to  the 
estimated  eeff. 

Generally,  for  L  one  or  two  wavelengths  long,  the  other  two 
values  will  be  very  far  off. 


B.  FINDING  THE  TRANSITION  PARAMETERS 


We  now  turn  to  the  rest  of  the  de-embedding  process. 

Using  identical  loads,  do  the  Reflect  measurements,  finding 
and  rMp.  It  is  not  necessary  to  know  the  load  value,  only  that 
it  is  the  same  for  the  two  reflect  measurements.  Taking  identi¬ 
cal  lengths  off  the  ends  of  the  TRL  transitions  provides  the 
simplest  of  loads. 


al 

transition 

b2 

al 

transition 

— ► 
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rL  rL  «- 

a0  bq 
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rMB 
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( b l/ a  1 )  A  “  rMA  ~ 


Ai(a2/t>2^  +  A2 
a3 (a2/b2^  +  ^ 


AirL  +  A2 
A3rL  +  1 


(14) 


define 


X1  -  RAl/RA3  ~  Al/A3 


and 


X2  “  RA2/RA4  ~  A2 


(15) 


(16) 


where  these  values  were  completely  determined  in  (7)  and  (9) . 
we  find  from  (14) 


Al= 


rMA  -  X2 


rL(1-rMA/X1) 


(17) 


where  all  values  are  now  known  except  rL 
Similarly, 


bl 

RBl/RB4  RB2/RB4 

• 

a2 

=  HB4* 

a! 

B 

RB3/RB4  1 

b2 

B 

1 

(bi/al) B  =  = 

rL 


B1  +  B2 (b2^a2 ) B 
b3  +  (b2^a2 ) B 


B1  +  B2  *  rMB 


B3  +  rMB 


(18) 


We  now  find  the  B  values,  hence  rL,  and  finally,  from  (17),  A-^. 


12 


From  measurements: 


Eliminating  the  unknown  rL  from  (17)  and  (18)  gives 

A1  <rMA  -  X2)^1  +  tB2/B1]-rMB) 

—  =  -  ,  all  known 

B1  (1  -  rMA/X1)-(B3  +  rMB) 


(22) 


Multiplying  (22)  by  (19)  yields  two  solutions  for 


where  a-^  is  a  positive  real  number. 


To  find  which  solution  is  to  be  used,  find  the  approximate  value 
of  A-^ •  If  an  ideal  open  circuit  were  used  for  the  reflect 
measurement,  rL  =  1.  If  the  open  circuit  was  at  an  angular  dis¬ 
tance  9^  before  the  reference  position,  Tt  =  c^201.  Substituting 
in  (17)  we  find  the  approximate  value  of  A-,. 

-  X2 

A1  «  -v,  . -  s  Alx  =  |alx|  *e^Tax,  all  known  (24) 

e  1)201  (1  _  rMA/Xl) 

Now  that  we've  calculated  the  phase  angle  that  A-j^  would  have 
with  an  ideal  open  circuit,  we  can  compare  that  phase  angle  with 
the  angles  derived  from  the  non-ideal  opens.  The  winner  is  the 
angle  closest  to  Tax.  This  is  shown  diagrammatical ly  in  Figure 
A 1  • 


Figure  Al.  Determining  the  correct  phase  of  A^ 
All  values  of  the  A  matrix  are  now  known. 
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C.  DE-EMBEDDING  THE  DEVICE 

The  device  is  inserted  at  the  reference  position,  that  is,  the 
junction  of  TRL  transitions  1  and  2  (Figure  A2) .  This  position 
need  not  be  at  the  center  of  the  line. 


Reference 

position 


TRL  Transition  A  TRL  Transition  B 

-line  length  Ld/2  -line  length  Ld/2 


Figure  A2.  Inserting  the  device  at  the  reference  position 


radb  can  be  written  to  include  the  negative  line  lengths 
J  -Ld/2  taken  by  the  DUT  from  the  TRL  transition  line 
lengths . 


The  R  matrix  for  negative  line  length  -Ld/2  is 


(25) 


so 

radb  ]  =  [  RA  ]  •  [  RN  ]  •  [  RD  ]  •  [  RN  ]  •  [  RB  ] 


~  RA4"RB4’  [  ARN  ]  •  | 

rd  ]  •  | 

rnb  ] 

(26) 

[-] 

[  “m  ]  • 

RA4 * RB4 

radb]  ' 

•  [*■»]  1 

(27) 
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The  values  RA4  and  RB4  have  not  yet  been  determined.  We  can 
write 


[  rab  ]  =  RA4*RB4  [  A  ] * [  B  ]  which  yields 

RAB4  =  RA4*RB4*(1  +  a3*B2)*  (28) 

Substituting  (28)  into  (27)  gives  the  result 


r  i 

(1  +  A3*B2)  r  1  -1 

r  I 

r  i -1 

L  r°  J 

=  p  *  L  ***  J • 

radb  • 

1  rnb  1 

(29) 

RAB4 

Examination  of  this  equation  shows  that  all  values  to  the  right 
of  the  =  sign  have  been  determined;  therefore  the  device  trans¬ 
mission  parameters  can  be  calculated. 


Finally,  the  4-port  device  R  parameters  can  be  converted  to  the 
device  S  parameters  by  the  standard  transformation  methods. 


_  _ 

r .  i  1 

R12  dr 

[_  SD  J  "  n 

•  dr  ~  Kll*k22  “  k12*k21 

R22 

1  “k21 

(30) 
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Appendix  B:  STRL4 ,  A  De-Embedding  Computer  Routine" 


The  computer  program  listed  in  the  following  pages  was  written  in 
TRUE  BASIC,  a  compiled  version  of  Basic  which  seems  well  suited 
for  mathematical  programming. 

The  first  part  of  the  program  is  specific  for  data  obtained  on  a 
Wiltron  360  network  analyzer.  It  calls  up  the  stored  S  parameter 
data  from  the  Transition,  added  Length,  Reflection,  and  DUT 
measurements  and  filters  out  the  garbage  associated  with  byte 
file  headings,  spaces,  etc.  The  four  Wiltron  data  files  are 
stored  with  file  names  T,L,R,and  any  name  for  the  device  under 
test.  All  four  files  are  automatically  given  the  extension  .DAT 
by  the  analyzer. 

Example:  For  TRL  de-embedding  if  we  denote  a  particular  measure¬ 

ment  set  "B,"  and  we  call  the  DUT  "fil,"  the  four  file  names 
would  be 
T_B.DAT 
LB . DAT 
R_B . DAT 
fil_B.DAT 

For  the  filtering  routine  to  work  with  the  Wiltron,  it  is  impor¬ 
tant  that  the  analyzer  is  set  to  display  linear  magnitude  and 
phase  information  such  as  "linear  polar",  and  that  the  markers 
are  off. 


If  data  is  obtained  from  another  analyzer  such  as  an  H.P.  8510,  a 
routine  should  be  written  so  that  data  for  N  different  frequen¬ 
cies  ends  up  as 

f (k) ,  k=l  to  N 

sll  (k,g)  ,  k=l  to  N,  g=l  to  4  (1-L,  2-+T,  3-R,  4->DUT) 

til (k,g) ,  "  "  ,  " 

•  •  •  • 

•  •  •  • 

t22 (k,g) ,  "  "  ,  »  " 


Copies  of  this  program  may  be  obtained  by  U.S.  companies  by 
writing  the  author.  Although  every  attempt  has  been  made  to 
assure  accuracy,  the  author  and  the  U.S.  government  do  not  assume 
any  responsibility  for  its  usage. 
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! STRL  (takes  4  Wiltron  byte  files  as  input  and  de-embeds  device 

S  parameters) 

input  prompt  "device  name,  device  length  (in.)=  ":dname$,  Ld 
input  prompt  "use  printer  (y  or  n)  ":print$ 
if  print$="y"  then 
open  #10:  printer 

print  #10:  dname  $  & " . TRL  (de-embedded) " , "  Ld=" ;Ld; "in. " 
print  #10: 
end  if 

input  prompt  "include  Ldb,  eeff  of  calib.  ?  (y  or  n)  " :incl$ 
option  nolet 
for  g=l  to  4 

if  g=l  then  open  #1:  name  "salb. wil" ,  access  input,  org  byte 

if  g=2  then  open  #2:  name  "sab. wil",  access  input,  org  byte 

if  g=3  then  open  #3:  name  "sga.wil" ,  access  input,  org  byte 

if  g=4  then  open  #4:  name  "dname$&" . wil" ,  access  input,  org  byte 

ask  #g:  filesize  lenl 
read  #g, bytes  lenl:file$ 

pl=pos [ f ile$,chr$ (10) &chr$ (10) &chr$ (49) ] -1  ! start  of  1st  freq. 

position 

p2=pos[file$,chr$(10) &chr$(12) ]  lend  of  data  position 

Ien2=p2-pl 

n=int (len2/137)  ! number  of  frequencies 
set  # g:  pointer  begin 

read  #g, bytes  pl:kill$  {removes  garbage  headings 

dim  buf$(141,10) ,c$(10) ,d$(10) ,f$(10) 
when  error  in 
for  k=l  to  n 

read  #g,  bytes  137:  a$  {reads  line  of  data 

b$=a$[2:137]  {removes  chr$(10)  symbol  before 

data 

for  j=l  to  10 

c$(l)=b$[01:03] 
c$ ( 2 ) =b$ [04:16] 
c$ (3 ) =b$ [17:31] 
c$ (4 ) =b$ [32:48] 
c$ (5) =b$ [49:63] 
c$ (6) =b$ [64:78] 
c$ (7) =b$ [79:93] 
c$ (8) =b$ [94:108] 
c$ (9) =b$ [109:124] 
c$ (10) =b$ [125:139] 

buf$(k, j)=c$(j) 
next  j 
next  k 
use 

print  using  "not  enough  data  on  freq  line  ###": k 
set  #g:  pointer  end 
end  when 

dim  f (401,4) , sll(401,4) , til (401, 4) ,sl2(401,4) ,cl2 (401,4) , 
s21 (401, 4) , t2 1 (401 , 4 ) ,  s22(401,4) ,t22(401,4) 
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for  k=l  to  n 
f (k,g)=val(buf$(k, 2)  ) 
sll (k,g) =val (buf$ (k, 3) ) 
til (k, g) =val (buf $ (k, 4) ) 
sl2 (k, g) =val (buf$ (k, 5) ) 
tl2 (k, g) =val (buf $ (k, 6) ) 
s21(k,g)=val(buf$(k, 7) ) 
t21 (k,g) =val (buf$ (k, 8)  ) 
s22 (k, g) =val (buf$ (k, 9)  ) 
t22 (k, g) =val (buf $ (k, 10) ) 
next  k 
close  #g 
next  g 

REM  g=l:Salb  g=2:Sab  g=3:Sga  g=4:Sadb 

REM:  example:  f(20,4)  is  the  20th  frequency  of  Sadb  data 
call  trl 


sub  TRL  !  TRL  calibration  and  measurement 
option  angle  degrees 
option  nolet 

La=.l!  added  line  (inches)  used  for  calibration 
re®  Ld=  length  (inches)  of  device  to  be  parameterized 
if  incl$=Hy"  then 

print  "F (GHz )  Ldb  eeff  Sll  All  S12  A12 

S21  A21  S22  A22" 

if  print$="y"  then  print  #10:  "F(GHz)  Ldb  eeff  Sll  All 

S12  A12  S21  A21  S22  A22" 


else 

print  "F (GHz)  Sll  All  S12  A12  S21  A21 

S22  A22" 

if  print$="y”  then  print  #10:  "F (GHz )  Sll  All  S12 

A12  S21  A21  S22  A22" 

end  if 

FOR  N=1  to  28 
FOR  G=1  TO  4 
fr=f (n,  1) 

!S  params  with  added  line 

!read  si , al , s2 , a2 , s3 , a3 , s4 , a4 

let  sl=sll(n,l) 

let  al=tll (n, 1) 

let  s2=sl2(n,l) 

let  a2=tl2(n,l) 

let  s3=s21(n,l) 

let  a3=t21 (n, 1) 

let  s4=s22(n,l) 

let  a4=t22(n,l) 
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call  StoT (sl,al,s2,a2 ,83,33,54,34, rlalb, tlalb, r2alb, t2alb, 
r3alb, t3alb, r4alb, t4alb)  !Ralb 
!S  params  without  added  line 
!read  Sl,al,s2,a2,s3,a3,s4,a4 
let  sl=sll(n,2) 
let  al=tll (n, 2 ) 
let  s2=sl2(n,2) 
let  a2=tl2 (n, 2) 
let  s3=s21(n,2) 
let  a3=t21(n,2) 
let  s4=s22(n,2) 
let  a4=t22(n,2) 

call  StoT(sl, al, s2 , a2 , s3 , a3 , s4 , a4 , rlab , tlab, r2ab, t2ab, r3ab, t3ab, 
r4ab,t4ab) !Rab 

! reflection  coeff.  of  transitions  and  B  with  identical  loading 

!read  ga, tga, gb, tgb 

let  ga=sll(n,3) 

let  tga=tll(n,3) 

let  gb=s22 (n, 3) 

let  tgb=t22(n,3) 

!S  parameters  with  device  inserted 

!read  sl,al,s2,a2,s3,a3,s4,a4 

let  sl=sll(n,4) 

let  al=tll (n, 4) 

let  s2=sl2(n,4) 

let  a2=tl2 (n, 4) 

let  s3=s21(n,4) 

let  a3=t21(n,4) 

let  s4=s22(n,4) 

let  a4=t22 (n, 4 ) 

call  StoT (si, al , s2 , a2 , s3 , a3 , s4 , a4 , rladb, tladb, r2adb, t2adb, 
r3adb, t3adb, r4adb, t4adb)  !Radb 

Ifind  matrix  [T]=[Ralb]*[Rab]-l 

call  inv(rlab, tlab, r2ab, t2ab, r3ab, t3ab, r4ab, t4ab, rilab, ailab, 
ri2ab, ai2ab, ri3ab, ai3ab, ri4ab, ai4ab) ! (Rab)-l 
call  mult (rlalb, tlalb, r2alb, t2alb, r3alb, t3alb, r4alb, t4alb, 
rilab, ailab, ri2ab, ai2ab, ri3ab,ai3ab,ri4ab,ai4ab, 
tl,atl,t2,at2,t3,at3, t4,at4) 

!find  line  coeff.  !  exp(-alpha*La) *exp(-j*beta*La)=gx*exp( j*tgx) 
!  first  solve  quad,  exp (-2*gamma*La) -exp (-gamma*La) * (T1+T4) 

+  (T1*T4-T2  *T3 ) =0 

call  add ( tl/2 , atl , t4/2 , at4 , b, angb) 

call  add ( tl*t4 , atl+at4 , -t2*t3 , at2+at3 , c , angc) 

call  add(l,0,-c/ (b*2) ,angc-2*angb,d,angd) 
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call  add(b,angb,b*sqr(d) , angb+angd/2 ,gl , tgl) 

call  add(b,angb, -b*sqr (d) , angb+angd/2, g2,tg2) 

if  tgl<=0  then  ! correct  solution  is  one  with  negative  exp. 

gx=gl 

tgx=tgl 

else 

gx=g2 
tgx=tg2 
end  if 

!gx*exp( jtg) =exp (-alpha * La) *exp(- j*beta*La) 
beta=-tgx/La 

eeff=(30*tgx/ (2*180*2. 54*La*fr)  )  ‘2 
Ldb=20*logl0 (Gx) /La 

! print  using  "f =##.###  Ldb=##.###  Eeff=## .###": fr,  Ldb,eeff 

! if  print$="y"  then  print  #10,  using  "f=##.###  Ldb =##.### 

Eef f=## . ##" :  fr,  Ldb,eeff 


! Find  Xl=Ral/Ra3,  X2=Ra2/Ra4 

call  add(  gx/t3,tgx-at3,-t4/t3,at4-at3,xl,txl) 
call  add ( 1/ (gx*t3 ) , -tgx-at3 , -t4/t3 , at4-at3 , x2 , tx2) 

a2=x2  lone  of  the  elements  of  [A]:  to  be  used  later 
ta2=tx2 

Iput  [Rab]  in  the  form  const* [abl , ab2 , ab3 , 1 ] 

call  ratios (r lab, t lab, r2ab,t2ab, r3ab, t3ab, r4ab, t4ab , abl , tabl , 
ab2 , tab2 , ab3 , tab3 ) 

lRb=R4b*[Bl,B2,B3,l] 

Isolve  for  B2/B1  =kl*exp ( j *tkl) 

call  add (ab2 , tab2 , -x2 , tx2 , ka , tka) 

call  add (abl , tabl , -x2*ab3 , tx2+tab3 , kb , tkb) 

kl=ka/kb 

tkl=tka-tkb 

Isolve  for  B3=b3*exp( j*tb3) 

call  add (ab3 , tab3 , -abl/xl , tabl-txl , ka , tka ) 

call  add ( 1 , 0, -ab2/xl , tab2-txl , kb, tkb) 

b3=ka/kb 

tb3=tka-tkb 

Isolve  for  Al*Bl=k3*exp( j*tk3) 

call  add (abl, tabl , -x2*ab3 , tx2+tab3 , ka , tka) 

k3=ka/kb 

tk3=tka-tkb 

Isolve  for  al*exp( j*tal) 

call  add(ga, tga, -x2 ,tx2 ,gl,tgl) 

call  add(  I,0,kl*gb,tkl+tgb,g2,tg2) 

call  add ( abl , tabl , -x2 *ab3 , tx2+tab3 , g3 , tg3 ) 

call  add(gb, tgb,b3,tb3,g4,tg4) 

call  add (1,0, -ga/xl , tga-txl , g5 , tg5 ) 

call  add (1,0, -ab2/xl , tab2-txl , g6 , tg6 ) 

al=sqr (gl*g2*g3/ (g4*g5*g6) ) 

taxl=(tgl+tg2+tg3-tg4-tg5-tg6)/2 
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tax2=taxl+180 


!find  if  tal=taxl  or  tal=tax2 

!find  approx  value  of  al,  given  that  load  was  an  exact  open 

tar=2*beta*lr/2  ! increased  load  angle  caused  by  missing  line 

call  add (ga , tga , -x2 , tx2 , anum, tanum) 

call  add (1,0, -ga/xl , tga-txl , aden, taden) 

ax=anum/aden 

tax=tanum-taden-tar 

if  abs (mod ( tax , 2*180) -mod (taxi , 2 *180) ) <180/2  then 
tal=taxl  else  tal=tax2 

a3=al/xl 

ta3=tal-txl 


a4  =  l 
ta4  =  0 

Isolve  for  Bl=bl*exp ( j *tbl) 

bl=k3/al 

tbl=tk3-tal 

! solve  for  B2=b2*exp ( j *tb2 ) 

b2=bl*kl 

tb2=tbl+tkl 


b4  =  l 
tb4  =  0 

! [Radb]  is  actually  [RaRnRdRnRb]  where  [Rn]  is  matrix  for  -Ld/2 
line  length 
! [Rn]  matrix  is: 
rnl=exp (alpha* Ld/2 ) 
tnl=beta*Ld/2 
rn2=0 
tn2=0 
rn3=0 
tn3  =  0 

rn4=exp ( -alpha* Ld/2 ) 
tn4=-beta*Ld/2 

!  find  [ Ran]  = [ Ra ] * [ Rn] =R4a* [A1,A2,A3,1]  and 
[ Rnb] = [Rn] * [Rb] =R4b* [ B1 , B2 , B3 ,1] 

!  [Rd]=[Ran] -1  * [Radb] * [Rnb] -1 

call  mult  (al ,  tal ,  a2  ,  ta2  ,  a3 ,  ta3 , 1, 0 ,  rnl ,  tnl ,  m2 ,  tn2  ,  m3  ,  tn3  , 
rn4 , tn4 , rani , tanl , ran2 , tan2 , ran3 , tan3 , ran4 , tan4 ) 
call  mult (rnl , tnl , m2 , tn2 , m3 , tn3 , rn4 , tn4 ,bl , tbl , b2 , tb2 , 

b3 , tb3 , b4 , tb4 , rnbl , tnbl , rnb2 , tnb2 , rnb3 , tnb3 , rnb4 , tnb4 ) 
i  find  [Rnb]-1  and  [Ran]-1 

call  inv (rnbl , tnbl , rnb2 , tnb2 , rnb3 , tnb3 , rnb4 , tnb4 , 

nbl , tinbl , inb2 , tinb2 , inb3 , tinb3 , inb4 , tinb4 ) 
call  inv ( rani , tanl , ran2 , tan2 , ran3 , tan3 , ran4 , tan4 , ianl , tianl , 
ian2 , tian2 , ian3 , tian3 , ian4 , tian4) 

! find  [Radb] *[Rnb-l]=Rbb 
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call  mult (rladb, tladb, r2adb, t2adb, r3adb, t3adb, r4adb, t4adb, 
inbl , tinbl, inb2 , tinb2 , inb3 , tinb3 , inb4 , tinb4 , 
rlbb, tlbb, r2bb, t2bb, r3bb, t3bb, r4bb, t4bb) 

!find  [Ran-1] *Rbb=Rdd 

call  mult ( ianl , tianl , ian2 , tian2 , ian3 , tian3 , ian4 , tian4 , 
rlbb, tlbb, r2bb, t2bb, r3bb, t3bb, r4bb, t4bb, 
rldd, tldd, r2dd, t2dd, r3dd , t3dd, r4dd , t4dd) 

! [Rd] =1/ (R4a*R4b)  *[Rdd],  where  1/ (R4a*R4b) = (A3*B3+1) /R4ab 

call  add(a3*b2/r4ab, ta3+tb2-t4ab, l/r4ab, -t4ab,de,tde) 

rdl=de*rldd 

rd2=de*r2dd 

rd3=de*r3dd 

rd4=de*r4dd 

tdl=tde+tldd 

td2=tde+t2dd 

td3=tde+t3dd 

td4=tde+t4dd 

!  Find  device  S  parameters  [Sd] 

call  TtoS (rdl , tdl , rd2 , td2 , rd3 , td3 , rd4 , td4 , sdll , tdll , sdl2 , tdl2 , 

sd2 1 , td2 1 , sd2 2 , td2  2 ) 

next  g 

if  incl$="y"  then 

print  using  "##.###  ##.##  #.##  ##.###  ####.# 

##.###  ####.#  ##•###  ####.#  ##.###  ####.#": 
f r , Ldb, eef f , sdll , tdll , sdl2 , tdl2 , sd21 , td21 , sd22 , td22 
if  print$="y"  then  print  #10,  using  »##.###  ##.##  #.##  ##.### 

####.#  ##.###  ####.#  ##.###  ####.#  ##.###  ####.#": 
f r , Ldb , eef  f , sdll , tdll , sdl2 , tdl2 , sd2 1 , td2 1 , sd2  2 , td2  2 

else 

print  using  ”##.###  ##.##  ##.##  ####.#  ##.###  ####.# 

##.###  ####.#  ##.###  ####.#••: 
f r , sdll , tdll , sdl2 , tdl2 , sd2 1 , td2 1 , sd22 , td22 
if  print$="y"  then  print  #10,  using  "##.###  ##.##  ##.##  ####.# 

##.###  ####•#  ##.###  ####.#  ##.###  ####.#": 

fr , sdll , tdll , sdl2 , tdl2 , sd21 , td21 , sd22 , td22 

end  if 
next  n 
end  sub 
end 

! *************  external  subroutines  follow  *************** 

SUB  MULT (al,tl,a2,t2,a3,t3,a4,t4,bl,ql,b2,q2,b3,q3,b4,q4, 

Cl,rl,c2,r2,c3,r3,c4,r4) 
call  add(al*bl, tl+ql,a2*b3 , t2+q3 ,cl, rl) 
call  add (al*b2 , tl+q2 , a2*b4 , t2+q4 ,c2 , r2 ) 
call  add (a3*bl , t3+ql , a4*b3 , t4+q3 , c3 , r3 ) 
call  add (a3*b2 , t3+q2 , a4*b4 , t4+q4 , c4 , r4 ) 
end  sub 

SUB  RTP(x,y,u,i)  !u*exp( j*i) =x+j (v2)  converts  x,y  to  u*expi 
option  angle  degrees 
u=sqr (x*2+y*2) 
if  y=0  then 
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if  x>=0  then  i=0  else  i=180 
end  if 

if  y<>0  then 

if  x=0  then  i=180/2*y/abs (y) 
if  x>0  then  i=atn(y/x) 
if  x<0  then  i=atn(y/x) +180 
end  if 
end  sub 

SUB  ADD (ul , il , u2 , i2 , u , i )  ! u*exp ( j *i) =ul*exp ( j *\1) +u2*exp ( j *i2 ) 

option  angle  degrees 

x=ul*cos ( il) +u2*cos ( i2 ) 

y=ul*sin(il)+u2*sin(i2) 

u=sgr (x*x+y*y) 

if  y=0  then 

if  x>=0  then  i=0  else  i=180 
end  if 

if  yoo  then 

if  x=0  then  i=180/2*y/abs (y) 
if  x>0  then  i=atn(y/x) 
if  x<0  then  i=atn(y/x)+180 
end  if 

!  i=360/ (2*180) *i 

end  sub 

SUB  INV (ul , il , u2 , i2 , u3 , i3 , u4 , i4 , vl , j 1 , v2 , j  2 , v3 , j  3 , v4 , j  4 ) 

rem  u's,  i's,  elements  of  matrix,  v's,j's  are  inverse 

call  add(ul*u4, il+i4,-u2*u3, i2+i3,d, id) 

vl=u4/d 

jl=i4-id 

v2=-u2/d 

j2=i2-id 

v3=-u3/d 

j3=i3-id 

v4=ul/d 

j4=il-id 

end  sub 

SUB  StoT(sl,al,s2,a2,s3,a3,s4,a4,tl,bl,t2,b2,t3,b3,t4,b4) 

!S  to  T  params 

call  add(sl*s4,al+a4,-s2*s3,a2+a3,d,td) I  find  delta  =d*exp(jtd) 

tl=-d/s3 

bl=td-a3 

t2=sl/s3 

b2=al-a3 

t3=-s4/s3 

b3=a4-a3 

t4=l/s3 

b4=-a3 

end  sub 

SUB  TTOS ( tl , al , t2 , a2 , t3 , a3 , t4 , a4 , si , bl , s2 , b2 , s3 , b3 , s4 , b4 ) 

!T  to  S  params 

call  add(tl*t4,al+a4,-t2*t3,a2+a3,d,td)  !delta=d*exp( jtd) 
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sl=t2/t4 

bl=a2-a4 

s2=d/t4 

b2=td-a4 

s3=l/t4 

b3=-a4 

s4=t3/t4 

b4=a3-a4+180 

end  sub 

SUB  Root (al , tal , b , tb , c , tc , xl , txl , x2 , tx2 ) 

Isolves  al (X~2) +B(X) +C=0  where  al=al*exp( jta) , . .etc 
! 2  roots  are  xl*exp( jtxl) ,  x2*exp(jtx2) 
if  al<>0  then 

call  add (1,0, -4*al*c/b" 2 , tal+tc-2*tb , e , te) 
call  add(l,0,  sqr(e) ,te/2, fl,tfl) 
call  add(l,0,-sqr(e) ,te/2,f2,tf2) 
xa=-b*fl/ (2*al) 


txa=tb-tal+tf 1 

xb=-b*f2/(2*al) 

txb=tb-ta 1+t  f 2 

if  abs ( xa ) >abs ( xb )  then 

xl=xa 

txl=txa 

x2=xb 

tx2=txb 

end  if 

if  abs ( xb ) >abs ( xa )  then 

xl=xb 

txl=txb 

x2=xa 

tx2=txa 

end  if 

end  if 

if  a 1=0  then 

xl=-c/b 

txl=tc-tb 

x2=xl 

tx2=txl 

end  if 

end  sub 

SUB  RATIOS ( rl , trl , r2 , tr2 , r3 , tr3 , r4 , tr4 , al , tal , a2 , ta2 , a3 , ta3 ) 

! forms  [R1,R2,R3,R4]  matrix  into  R4* [Al, A2 , A3 , 1]  matrix 

al=rl/r4 

tal=trl-tr4 

a2=r2/r4 

ta2=tr2-tr4 

a3=r3/r4 

ta3=tr3-tr4 

end  sub 
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